Spin dynamics and spin freezing 
in the triangular lattice antiferromagnets FeGa2S4 and NiGa2S4 
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Magnetic susceptibility and muon spin relaxation (^SR) experiments have been carried out on 
the quasi-2D triangular- lattice spin S = 2 antiferromagnet FeGa2S4. The /iSR data indicate a 
sharp onset of a frozen or nearly-frozen spin state at T* = 31(2) K, twice the spin-glass-like freezing 
temperature Tf = 16(1) K. The susceptibility becomes field dependent below T* , but no sharp 
anomaly is observed in any bulk property. A similar transition is observed in /iSR data from the 
spin-1 isomorph NiGa2S4. In both compounds the dynamic muon spin relaxation rate \d{T) above 
T* agrees well with a calculation of spin-lattice relaxation by Chubukov, Sachdev, and Senthil in 
the renormalized classical regime of a 2D frustrated quantum antiferromagnet. There is no firm 
evidence for other mechanisms. At low temperatures Xd{T) becomes temperature independent 
in both compounds, indicating persistence of spin dynamics. Scaling of \d{T) between the two 
compounds is observed from ~T/ to ~1.5r*. Although the /xSR data by themselves cannot exclude 
a truly static spin component below T* , together with the susceptibility data they are consistent 
with a slowly-fluctuating "spin gel" regime between T/ and T*. Such a regime and the absence of 
a divergence in A£j(r) at T* are features of two unconventional mechanisms: (1) binding/unbinding 
of Z2 vortex excitations, and (2) impurity spins in a nonmagnetic spin-nematic ground state. The 
absence of a sharp anomaly or history dependence at T* in the susceptibility of FeGa2S4, and the 
weakness of such phenomena in NiGa2S4, strongly suggest transitions to low-temperature phases 
with unconventional dynamics. 

PACS numbers: 75.10.-b, 75.10.Hk, 75.10.Jm, 76.75. 



I. INTRODUCTION 



Among geometrically frustrated magnets the simplest 
example is the two-dimensional (2D) triangular-lattice 
Heisenberg antiferromagnet (THAFM) with isotropic 
spin-spin interactions. This is the system for which 
a spin liquid state in more than one dimension was 
first proposed by Anderson^ as a quantum disordered 
state in which long-range magnetic order is destroyed. 
It is now beheved that the spin- 1/2 2D THAFM with 
nearest-neighbor coupling orders at T = with a 120° 
spin structure. A number of treatments of the 2D 
THAFM at nonzero temperatures have been reported.'^ 
Phase transitions associated with the value of the spin S, 
topological defects, interactions beyond nearest neigh- 
bors, additional terms in the Hamiltonian, etc., have been 
studied intensively. One such transition involves binding 
and unbinding of topological excitations, the so-called 
Z2 vortices. Integer spins and additional biquadratic 
nearest-neighbor coupling are predicted to yield spin- 
nematic phases with no dipolar ordering. ^^'^^ Coupling 
beyond nearest neighbors may lead to a quantum spin 
disordered ground state^"^ or a state with broken C3 
symmetry.^^ds 1^ j^^^g proved difhcult to find experimen- 
tal evidence for many of these transitions, although Z2 
vortex binding has been proposed for a number of candi- 



date 2D THAFM systems^^ and a recent theory^^ invokes 
impurity-spin interactions in a nonmagnetic spin-nematic 
low-temperature phase. 

The quasi-2D triangular-lattice antiferromagnetic in- 
sulator NiGa2S4 was characterized by Nakatsuji and co- 
workers^*"'^® and found to possess a number of unusual 
properties that have attracted considerable attention, 
both experimentallyi^"^'* and theoretically ^^^^'^^^^'^^^^i 
It is a 2D THAFM with spin 5 = 1 and very httle 
spin anisotropy. Third-nearest-neighbor antiferromag- 
netic exchange is dominant, with a weak nearest- neighbor 
ferromagnetic interaction^^'^^ as suggested by calculated 
superexchange pathways.*^ A phase transition, possibly 
influenced by impurity effects, is indicated by a cusp 
and weak bifurcation between field-cooled (FC) and zero- 
field-cooled (ZFC) dc magnetic susceptibilities x<ic{T) at 
T* = 8.5-9 

A drastic slowing of Ni^+ spin fluctuations as T* is 
approached from above is observed in NMR,^° ESR,^^'^^ 



and muon spin relaxation (/xSR) 



23,25-28 



experiments. 



with a frozen or nearly-frozen Ni^"'' spin state below T* 
characterized by an order-parameter-like temperature de- 
pendence of the (nearly) frozen moment and strong spin 
fluctuations down to 25 mK. The frequency-dependent 
ac magnetic susceptibility Xa.c{T) suggests a spin-glass- 
like transition at a lower temperature Tf = 2.2-2.7 K.^^ 
Based on these results, a viscous "spin-gel" or extended 
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thermodynamic critical regime has been proposed be- 
tween ~Tj and T*}^'"^^ A similar broad fluctuating 
regime is observed in NaCr02.^^ 

Perhaps the most mysterious property of NiGa2S4 is 
the magnetic specific heat Cm, which is independent of 
applied magnetic field up to 7 T and exhibits a tem- 
perature dependence at low temperatures.^'' This implies 
linearly-dispersing 2D excitations that do not couple to 
the field, which seems hard to reconcile with the observa- 
tion that the low-temperature muon spin relaxation rate 
is rapidly suppressed by applied field. 

The isostructural insulator FeGa2S4,"'^*'"'^^ in which 
Fe^"*" is in the ilg^g high-spin S = 2 configuration, is 
also a candidate for a 2D THAFM. Both FeGa2S4 
and NiGa2S4 exhibit FC-ZFC bifurcation at ^\ew\/10, 
where 6w is the (negative) Weiss temperature from the 
paramagnetic-state susceptibility at high temperatures. 
Normally bifurcation occurs at a spin freezing tempera- 
ture Tf. Thus the Ramirez parameter^"* 9w/Tf is large 
(> 10), indicating strong frustration. In both compounds 
Cm exhibits an unusual two-peak temperature depen- 
dence with no sign of a phase transition. As in NiGa2S4, 
Cm{T) in FeGa2S4 follows a field-independent law at 
low temperatures. Some properties of the two materials 
differ significantly, however. In NiGa2S4 Xac(7') exhibits 
a frequency dependence below Tj,^^ whereas in FeGa2S4 
bifurcation, much stronger than that at T* in NiGa2S4, 
is observed at Tf « 16 K,^^ indicating a spin-glass-like 
transition at this temperature. In the Nii^a;FexGa2S4 
series the specific heats scale with 9w in the region 
but not at higher temperatures.^* 

To date FcGa2S4 has been studied using only a few mi- 
croscopic techniques. In Nii_j;Fej;Ga2S4, 0.01 < x < \, 
^'^Fe Mossbauer-effect (ME) measm-cmcnts^°'^^ indicate 
spin freezing at a temperature T* that varies smoothly 
between 12(1) K for x = 0.01 and 33(1) K « 2T/ for 
X = 1,^^ with an order-parameter- like temperature de- 
pendence of the frozen moment and strong spin fluctu- 
ations below T* similar to the /xSR results in NiGa2S4. 
Initial /iSR experiments in FcGa2S4'*''' revealed a mag- 
netic transition at and strong muon relaxation below 
T* w 30 K, confirming the ME results. It was noted 
that the relation of T* to Cm (T) is not the same in the 
two compounds: in NiGa2S4 T* lies just below the low- 
temperature peak in Cm{T), whereas in FcGa2S4 T* is 
found between the two peaks. As a result Cm(T) does 
not scale with T* , and the relation between spin freezing 
and entropy release is unclear. 

This paper (1) reports further measurements of Xdc(I') 
and results of /zSR experiments in FeGa2S4, which were 
undertaken to clarify the nature of spin freezing and to 
probe Fe^+ spin dynamics, and (2) compares /iSR and 
other data from FeGa2S4 with those from NiGa2S4. The 
susceptibility data confirm earlier results, -^'^ and extend 
measurements to applied fields in the range 0.1-2 T. As 
previously noted, there is no anomaly in the muon re- 
laxation rate at Tf and also no sharp anomaly in Xdc 
at T*; the latter is an extremely unusual feature of this 



compound,™ not shared by NiGa2S4. 

The /xSR results from FeGa2S4 provide new insight 
into spin freezing and dynamical fiuctuations in these 
compoimds. The data reveal an abrupt transition at 
T* « 31 K to a phase with frozen or nearly- frozen spins 
and strong spin fluctuations but with almost no signature 
in bulk measurements. The agreement with the ME re- 
sults rules out perturbation by the muon electric charge*^ 
as the origin of this behavior. 

The theory of spin-wave excitations in a 2D frus- 
trated quantum antiferromagnet has been treated by 
Chubukov, Sachdev, and Senthil^s-so (hereafter CSS) 
and others. ^^'^^ Comparison of muon relaxation rates 
above T* in both compounds with the CSS result for the 
spin-lattice relaxation rate of a probe spin (nucleus or 
muon) yields values of the id exchange interaction in re- 
markably good agreement with previous studies. ^^'^^'^^'^^ 
There is no firm evidence for muon relaxation due to 
other mechanisms. ^'^'^'^^ 

The muon relaxation rates in FcGa2S4 scale with those 
in NiGa2S4^^ over a wide temperature range from ~T/ 
to ~1.5T*. Although the /zSR data by themselves do 
not rule out a truly static spin component, as found 
in conventional magnetic phases, the spin dynamics re- 
vealed by juSR together with the absence of a suscep- 
tibility anomaly at T* are consistent with an extended 
critical regime or slowly-fluctuating "spin-gel" state in 
both compounds below T* . 



II. EXPERIMENTAL PROCEDURE AND 
RESULTS 

A. Magnetic susceptibility 

A powder sample of FeGa2S4 was prepared as de- 
scribed previously. ^^'^^'^* The bulk dc susceptibility was 
measured using a Quantum Design Magnetic Properties 
Measurements System over the temperature range 2- 
300 K for applied magnetic flelds in the range 0.1-6 T. 
For FC measurements the magnetic fleld was set and the 
sample was cooled to 2 K before collecting the data, 
while for ZFC measurements the sample was cooled in 
zero field, the field was set at 2 K, and data were taken 
upon warming. After each measurement the sample was 
warmed to 350 K in zero field to quench any magnetic 
order. 

The results are shown in Fig. 1. Consistent with a pre- 
vious report, for applied field jioH = 0.1 T bifurcation 
between the FC and ZFC data is observed below 16(1) K, 
suggesting spin-glass-like freezing at this temperature. 
The data indicate a strong suppression of the spin freez- 
ing temperature with applied field, similar to that ob- 
served in Fe-rich (Ni,Fe)Ga2S4 alloys,^''' to about 10 K 
at 2 T; this is also consistent with spin- glass-like freez- 
ing. As in NiGa2S4,^^ the susceptibility is suppressed by 
field (i.e., a negative nonlinear susceptibility sets in) for 
temperatures ^ T* obtained from ME and /xSR experi- 
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FIG. 1. (Color online) Temperature dependence of the dc 
magnetic susceptibility Xdc of FeGa2S4 for various applied 
magnetic fields. Filled symbols: field-cooled (FC) data. Open 
symbols: zero-field-cooled (ZFC) data. T/: low-field spin 
freezing temperature. T*: onset temperature of quasistatic 
muon spin relaxation [cf. Fig. 3(a)]. 



ments [cf. Fig. 3(a) below], but there is no sharp feature 
in Xdc{T) at T*. As noted in Sec. I, the bifurcation in 
FeGa2S4 is much stronger than that in NiGa2S4. 



B. Muon spin relaxation 

The /iSR technique is a sensitive probe of static and 
dynamic magnetism in magnetic materials. ^^'^^ In fiSR 
experiments spin-polarized positive muons are implanted 
in the sample and come to rest at interstitial sites. Each 
muon precesses in the sum of the local field Bioc(i) due 
to its magnetic environment and any applied field, and 
decays with a mean lifetime = 2.197 fis into a positron 
and two neutrinos. The positron is emitted preferentially 
in the direction of the muon spin at the time of decay, 
so that detection of decay positrons permits determina- 
tion of the evolution of the muon spin polarization. The 
resulting signal can be observed up to times of several 
Tfj. (typically 10-15 /is) and is conceptually similar to the 
free induction decay of NMR,^^ although the detection 
technique is very different. 

The onset of a quasistatic^^ component (Bioc) of 
Bioc(i) is expected at a magnetic phase transition. In 
general muon spin relaxation [decay of the ensemble 
muon spin polarization P{t)] is due to a combination of 
two relaxation mechanisms:^^'^^ 

• quasistatic relaxation, due to dephasing of muon 
spin precession in (Bioc) if its magnitude is inho- 
mogeneously distributed, and 

• dynamic (spin-lattice) relaxation of the muon spin 
due to thermal fluctuations of Bioc(i). 

These are of course the inhomogeneous and homogeneous 
relaxation mechanisms of NMR. 



It is straightforward to separate the relaxation rates 
associated with these processes in /zSR experiments in 
zero and weak longitudinal magnetic field iJ^,^^'^^ pro- 
vided that the dynamic muon relaxation is slower than 
the quasistatic muon relaxation. Then P{t) has a char- 
acteristic two-component form that can be modeled by 

P{t) = {l-ru)P,{t)+VdPd{t). (1) 

Here the first and second terms on the right-hand side de- 
scribe quasistatic and dynamic relaxation processes that 
control the evolution of P{t) at early and late times, 
respectively,^'''^'' and rjd is the fraction of the initial muon 
spin polarization that is relaxed dynamically, i.e., that 
survives at late times after Pq{t) 0. Observation of 
such two-component behavior indicates the presence of 
quasistatic magnetism, with or without long-range or- 
der. In a randomly-oriented powder sample % = 1/3 for 
(Sloe) > fJ-oHL and ?7d -> 1 for (Sioc) < PI'qHl,^° so that 
for zero or low Hl a sudden decrease of rjd with decreas- 
ing temperature signals the onset of a quasistatic local 
field due to a magnetic phase transition. 

/iSR experiments were carried out on a powder sample 
of FeGa2S4 using the M20 beam line at TRIUMF, Van- 
couver, Canada. The asymmetry A(t) in positron count 
rate, which is proportional to P(t), was measured using 
the standard time-differential ^SR technique. ^^'^^ Data 
were taken from 2 K to 260 K in a weak longitudinal 
field hoHl ~ 2 mT to decouple nuclear dipolar fields,^" 
and also for ^qHl between 2 and 100 mT at 1.7 K. 

Representative late-time muon asymmetry data for 
IIqHl ~ 2.02 mT are shown in Fig. 2(a). As previ- 
ously reported,**^ the initial asymmetry decreases rapidly 
between 40 and 20 K. This indicates^" the onset of 
quasistatic Fe^^ spin freezing (reduction of ry^) as dis- 
cussed above. As noted in Sec. I, ME studies^°''^^ of 
Nii_^Fea;Ga2S4 give evidence for quasistatic Fe^+ spins 
below T*, consistent with this loss of asymmetry. Each 
technique places a lower limit on the correlation time Tc 
of (Bloc). In ME and /iSR experiments these limits are 
-10"^ s and -lO-'^-lO'^ s (cf. Sec. IIIB), respectively, 
so that the /iSR results are a somewhat stronger. 

Unfortunately the early-time quasistatic muon relax- 
ation could not be observed at low temperatures in 
FeGa2S4,'''^ because the initial decay time was shorter 
than the spectrometer "dead time" between a muon 
stop and the earliest detection of the decay positron 
(^10 ns). Thus (i?ioc) is large and broadly distributed. 
The spectrometer dead time puts a lower limit of 
roughly 200 yits"^ on the quasistatic muon relaxation 
rate Aq, corresponding to a quasistatic local field dis- 
tribution width (i?ioc)rms = ^q/l^i ^ 0.2 T; here 7^ = 
851.56 /is~^ T~^ is the muon gyromagnetic ratio. The 
average local field (i3ioc)av may be considerably larger 
than this; in NiGa2S4, where the early-time asymmetry 
was observable, (i?ioc)av ^ 5 (Sioc)rms (Refs. 23 and 25). 

The asymmetry data were fit to a relaxation function 
of the form 

Ait) = AoP{t) , (2) 
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FIG. 2. (Color online) Representative late-time muon asym- 
metry data (positron count asymmetry vs time) in FeGa2S4. 
(a) Asymmetry vs time at various temperatures for longitudi- 
nal field hoHl = 2.02 mT. (b) Asymmetry vs time for various 
values of hoHl at temperature T — 1.7 K. Curves: fits to 
Eq. (2). 

where Aq is the initial count-rate asymmetry and P(t) is 
given by Eq. (1) with 

Pdit) ^ exp[-{Xdtf] , (3) 

the stretched-exponential form with relaxation rate 
and stretching power /? < 1. This is a convenient 
parametrization of sub-exponential relaxation due to an 
inhomogeneous distribution of locally exponential relax- 
ation processes. ^'^'^^ is a characteristic relaxation rate 
(not the average), and /3 controls the width of the dis- 
tribution, which becomes larger for smaller /3. Represen- 
tative fits are shown in Fig. 2. The stretched-exponential 
form was also used to analyze late-time data from /iSR 
experiments in NiGa2S4 using powder samples^^'^^'^^'^* 
and a mosaic of oriented single crystals with nor- 
mal to the Ni planes. The origin of this relaxation- 
rate inhomogeneity is not clear, since structural studies 
of these materials found nearly perfect triangular NiS2 
lattice layers with no distortion. It does not seem to 
be the stacking faults inferred from NQR,^°'^* since the 
bimodal distribution this would produce is not observed. 
We note that spin disorder and a short spin-spin correla- 
tion length are found in neutron scattering experiments 
on NiGa2S4."'^^ 



FIG. 3. (Color online) (a) Temperature dependence of dy- 
namic muon spin relaxation rate Ad in weak longitudinal ap- 
plied field Hl in FeGa2S4 and NiGa2S4. Circles: FeGa2S4 
powder, hqHl = 2.02 mT. Squares: NiGa2S4 powder, 
^oHl ~ 2.0 mT. Triangles: NiGa2S4, mosaic of oriented 
single crystals, || c, fioH^ = 2.0 mT. r/(FeGa2S4): freez- 
ing temperature from FC-ZFC bifurcation in Xdc(r) (Fig. 1). 
r/(NiGa2S4): freezing temperature from onset of frequency- 
dependent Xac(r) (Refs. 18 and 29). T*: transition temper- 
atures from /iSR data. Inset: FeGa2S4, late-time asymmetry 
fraction rjd [Eq. (1)] (filled circles) and stretching power f5 
[Eq. (3)] (open circles), (b) Normalized muon relaxation 
rate Xd{T)/Xd{T*) vs normalized temperature T/T*. Sym- 
bols as in (a). 



The temperature dependencies of A^; for fioH^ ~ 2 mT 
in FeGa2S4 and NiGa2S4 are shown in Fig. 3(a). ^'^ Sharp 
cusps in Xd{T) are observed at T* = 31(2) K and 
8.5(1.0) K in FeGa2S4 and NiGa2S4, respectively. In both 
compounds Xd{T) increases by more than two orders of 
magnitude over a wide temperature range as T — > T* 
from above. This is clear evidence for critical slowing 
down of magnetic fluctuations: as in NMR, A^ is propor- 
tional to the spin correlation time in the motional nar- 
rowing limit appropriate to the paramagnetic state. 
From Fig. 3(a), in FeGa2S4 T* is about twice the spin 
freezing temperature Tf. 

For T < T* Xd{T) decreases with decreasing temper- 
ature in both compounds but saturates and shows no 
sign of vanishing as T ^ (down to T/T* 0.003 in 
FeGa2S4*^). Such persistence of strong spin dynamics to 
low temperatures is a general feature of geometrically 
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frustrated magnets^^ and is not well understood, but 
seems to indicate a zero-energy singularity in the den- 
sity of excited states. In FeGa2S4 there is no anomaly 
in Xd(T) at Tf. This is very surprising, and is not under- 
stood: the observed strong FC-ZFC bifurcation is usu- 
ally considered good evidence for spin freezing, which in 
turn would be expected to reduce the relaxation rate. In 
NiGa2S4 a minimum in Xd{T) is seen near Tf, but data 
above and below this temperature were taken in different 
cryostats and an experimental artifact cannot be ruled 
out. 

The temperature dependencies of the late-time asym- 
metry fraction r]^ and the stretching power (3 in FeGa2S4 
are shown in the inset to Fig. 3(a). At low tempera- 
tures rjd « 0.3-0.5, close to the value 1/3 expected for 
HqHl ^ {B\oc}, and then rises rapidly to ~1 at T* . This 
is strong evidence for frozen or nearly-frozen Fe^"*" spins 
below T* , since it is the expected behavior if (Bioc) sets 
in suddenly at this temperature. 

Below T* /3 is considerably smaller than 1, although it 
is not accurately determined by the data; it was necessary 
to constrain /3 in this region to obtain consistent fits.'*^ 
We set a minimum /3 = 0.4, but the behavior of Ad(T) is 
not very sensitive to the specific choice. For T > T* f] 
tends to be larger but still < 1. The stretching has been 
found to be less pronounced above T* at other applied 
fields in both FeGa2S4 and NiGa2S4.23'^^'^^ 

Figure 3(b) shows the dependence of the normalized 
muon relaxation rate \d{T)/\d{T*) on the normalized 
temperature T/T* in FeGa2S4 and NiGa2S4. It can be 
seen that the data scale from to ^1.5T*. This be- 
havior and the possibility of slow quasistatic fluctuations 
below T* (cf. Sec. IIIB) suggest an extended critical re- 
gion, followed by disordered freezing at Tf}^'"^^ Scahng 
does not hold for the transition temperatures themselves, 
however: the ratio T*/\ew\ is 0.19(1) for FeGa2S4 and 
0.11(1) for NiGa2S4. 

The field dependencies of Xd and rjd in NiGa2S4 and 
FeGa2S4 for T ^ 2 K are shown in Fig. 4. It can 
be seen that in both compounds the muon relaxation 
is suppressed by field. For FeGa2S4 varies roughly 
as iJi"^'^, and is suppressed by an order of magnitude 
for iiqHl ~ 10 mT. For fields an order of magnitude 
higher than this rjd remains in the neighborhood of 0.4- 
0.5, close to the value 1/3 expected for jioHl <C (i3ioc)-^° 
In NiGa2S4 the field dependence of Ad is slower but still 
appreciable: Xd^HoHL) « 0.1Ad(~l mT) at jioHl ~ 
100 mT, where the increase of r/^ is just beginning. This 
demonstrates that for both compounds the suppression 
of Xd is not due to decoupling by , which is significant 
only as rjd ~^ 1.^° It is also not due to the glassy dynam- 
ics mechanism of Keren et al.,^^ which is not applicable 
when fioHL ^ Bioc. The suppression occurs for fields 
two to three orders of magnitude smaller than ksT* / jiB, 
which is far too small to reduce the 3(i-spin fluctuation 
amplitude by simple alignment of the id spins. A shift 
of spectral noise power to higher frequencies is required, 
which must be due to some less direct mechanism. The 
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FIG. 4. (Color online) Dependence of dynamic muon spin re- 
laxation parameters on longitudinal field Hl in FeGa2S4 and 
NiGa2S4 at low temperatures (T ^ T*). Circles: FeGa2S4 
powder, T = 1.7 K. Squares: NiGa2S4 powder, T = 2.1 K. 
Triangles: NiGa2S4 oriented single- crystal mosaic, T = 2.3 K. 
(a) Relaxation rate A^. (b) Late-time asymmetry fraction r]d- 



nature of this mechanism, and why it results in a power- 
law field dependence for FeGa2S4, are unknown. 

In NiGa2S4 rjd increases from its low- field value for 



> 



cjp(0)/7^ ~ 250 mT, where ^^(0) is the ob- 



served spontaneous muon precession frequency at T = 
Q 23,25 ^j-^ increase is expected at precisely this field, as 
Hl decouples the muon spin polarization from the qua- 
sistatic field. This quantitative agreement is strong ev- 
idence that the two-component structure of the muon 
relaxation function is not associated with separation of 
magnetically distinct phases. In that case rjd would be 
the relative fraction of the slowly-relaxing phase; it and 
the initial 100% muon spin polarization would both be 
field-independent. The stretched-exponential form of the 
late-time relaxation function indicates inhomogeneity of 
the muon relaxation rate, but there is no evidence for 
magnetic phase separation. 

Below 100 mT rjd varies somewhat more with field in 
FeGa2S4 than in NiGa2S4. This behavior is not under- 
stood, but might involve the significantly broader local 
field distribution in FeGa2S4. 
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III. DISCUSSION 
A. T > T*: comparison with the CSS result 

We first consider muon spin relaxation for T > T* .^"^ 
The CSS calculation of the probe-spin dynamic relax- 
ation rate 1/Ti due to spin fluctuations in a 2D frustrated 
quantum antiferromagnet, discussed in more detail in the 
Appendix, yields 

VT,.Cces(^)'A(J)%xp(T./T, ,4) 

in the renormalized-classical (RC) regime T <^ STrp^, 
with 

To = Airp, . (5) 

Here NqAq is a renormalized hyperfine coupling 
constant, Ps is the spin stiffness constant, and Cess 
is a numerical constant that can be estimated from 
Eqs. (4.29) and (4.7) in the Appendix. It is equal 
to 114.59 if one assumes Eq. (4.29) to be an equality, 
but this is an order-of-magnitude estimate. The the- 
ory assumes nearest-neighbor AFM coupling, so that 
for NiGa2S4, where third-nearest-neighbor interactions 
dominate,^ ^'^^ the parameters obtained are effective 
values.^® 

In general p|| and pj_ of ps for spin twists parallel 
and perpendicular, respectively, to the 2D plane are 
not equal. In the present experiments the samples 
were randomly-oriented powders except for the mosaic of 
NiGa2S4 single crystals, for which Ad(T) is not very dif- 
ferent from the powder. We therefore take Ps in Eq. (5) 
to be the orientation average (2p|| -I- p±)/3. 

We compare Eq. (4) to the muon relaxation rate A^, 
recognizing that the uncertain absolute value of l/Ti 
from Eq. (4) and the inhomogeneity make full agreement 
unlikely. Nevertheless the predicted T^exp(To/T) tem- 
perature dependence might be expected. Figure 5 is a 
semi- log plot of Xd/T^ vs 1/T for both compounds at 
temperatures above T*. For each compound the data 
extend down to the transition, i.e., the upper right-hand 
data points are for T w T* . 

For FeGa2S4 the straight line is a fit of Eq. (4) to the 
lowest temperature points. For NiGa2S4 a straight line 
fits all the data for both powder and mosaic samples. The 
resulting values of To yield l/27rp, = 2/To = 0.0068 K'^ 
and 0.0154 for FeGa2S4 and NiGa2S4, respectively 
(arrows labeled by symbols in Fig. 5), thereby justifying 
a posteriori the procedure of obtaining ps from Tq.^^ 

For NiGa2S4 all the data satisfy this condition, 
whereas for FeGa2S4 the data fall below the straight line 
at small 1/T. This dropoff might be due to the onset of 
thermally-activated muon diffusion. It should be noted, 
however, that it occurs for T « 27rps, where the RC 
calculation ceases to be valid, and that similar behavior 
was observed in "^Li NMR in the THAFM LiyRuOe in the 
temperature range 20-40 K,^^ where Li diffusion would 
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FIG. 5. (Color online) Dependence of Xd/T^ on inverse tem- 
perature 1/T in FeGa2S4 and NiGa2S4, T > T* . Straight 
lines: fits of Eq. (4). Labeled arrows: values of l/2n ps (see 
text). Curve: fit of the Z2 vortex theory^'^'' (Sec. IIIC) to 
data for FeGa2S4. 

not be expected. The situation is discussed further be- 
low. The curve is a fit of the Z2 vortex theory^'^'^ to the 
data, as described in Sec. IIIC. 

The spin stiffness constants are expected to be propor- 
tional to the exchange constant J; for a square lattice 
Ps = JS^ in the classical limit.^^ For the 2D THAFM py 
and p± have been calculated in the spin-wave approxi- 
mation, which yields^°'^^ 

Ps ^ (2/3)p|| + il/3)pi_ ^ 1 ~ 0.399/25 
J52 J52 y/3 ^ ' 

to first order in 1/S. Equation (6) and Ps obtained from 
the /iSR data give a "spin-wave" value Jsw of J, which 
can be compared with the value Jg^^ obtained^*^ from 
\6w\ = zJg„S{S + l)/3 for AF exchange between z = 6 
near neighbors. 

Experimental values of these and other spin-system 
parameters and are given in Table I together with the- 
oretical values of Ps/ J ■ It can be seen that for both 
compounds the values of Jsw and J^^, are in good agree- 
ment. This is strong evidence that their spin dynamics 
above T* are dominated by the spin-wave fluctuations 
treated by CSS. For NiGa2S4 there is also reasonable 
agreement with the values J3 = 32(7) K and 21 K of the 
dominant third-nearest-neighbor exchange constant ob- 
tained, respectively, from neutron scattering'^^ and the 
field dependence of the ESR resonance frequency well 
below T*.2i The spin stiffness constants for the two com- 
pounds are in rough and good agreement with the values 
of r* for FeGa2S4 and NiGa2S4, respectively, as expected 
for Halperin-Saslow modes^^ in spin-frozen states. 

The renormalized muon-Ni2+ hyperfine coupling con- 
stant NoAo/h estimated from Eq. (4) is expected to be 
of the order of the spontaneous T=0 muon precession 
frequency w^(0) = 7/j(i?ioc)av(0), since the same inter- 
action is involved in both cases. In FeGa2S4 w^(0) has 
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TABLE I. Experimental and theoretical values of spin-system parameters in NiGa2S4 and FeGa2S4. Tf. freezing temperature 
from Xac (NiGa2S4) and FC/ZFC bifurcation (FeGa2S4). T*: transition temperature from /xSR and ME data. 6w'- (negative) 
Weiss temperature from x{T)- i^m(O): spontaneous T=0 muon precession frequency (NiGa2S4 only). NoAq/H, Tq: estimated 

renormalizcd muon hypcrfinc coupling constant (Rof. 49) and characteristic temperature from fits of Eq. (4) to \d{T>T*) 
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not been measured (cf. Sec. 11 B), but in NiGa2S4 it is 
considerably larger than NqAq/Ti obtained from Eq. (4) 
(Table I). 

There are a number of uncertainties in the comparison, 
including the prcfactor Cess in Eq. (4) and the inhomo- 
geneity in the relaxation. Perhaps the least certain quan- 
tity in the CSS expression for 1/T\ is the value of ^/c,^® 
where ^ is the spin correlation length and c is the spin- 
wave velocity. Large discrepancies between calculated 
and observed correlation lengths have been reported. 
The situation is discussed in the Appendix. We conclude 
that the temperature dependence of in NiGa2S4 is in 
good agreement with the CSS result, but the; magnitude 
of the muon relaxation rate is not well understood. 

As shown in the Appendix, in the CSS treatment^^ 
the spin-lattice relaxation rate falls below the RC tem- 
perature dependence in the crossover region between the 
RC (T < 2-Kps) and quantum-critical (QC) (T » 2-Kps) 
regimes. Thus the high-temperature dropoff in the re- 
laxation rate in FeGa2S4 (Fig. 5) might be due to this 
crossover. However, the QC regime, for which the condi- 
tion Ps ^ is necessary, is suppressed for large spin [cf. 
Eq. (6)] and probably does not exist in FeGa2S4 [S = 2) 
because ps ~ 2J (Table I). Nevertheless, rapid non- 
universal short-range spin fluctuations may dominate at 
high temperatures^® and reduce Tc- The QC region might 
exist in NiGa2S4 [S = 1), where ps w 0.5J (Table I), 
but here there is no sign of a dropoff for T < 30 K 
{T/^-Kps ^ 0.46); apparently the crossover has not been 
reached. Above ~30 K Ad becomes too small to be mea- 
sured reliably. 

The ESR linewidth in NiGa2S4 above T* has been an- 
alyzed in terms of lifetime broadening by critical spin 
fluctuations. ^^'^^ The linewidth depends much less on 
temperature than the muon relaxation rate, and exhibits 
an anomaly at 23 K ~ 3T* that was interpreted as a 
crossover between critical regimes. It should be noted 
that the ESR experiments were carried out in a field 
of 20 T, corresponding to a Zeeman energy ~ 2kBT*; 
such a large field is expected to perturb the spin fluc- 
tuation spectrum in theoretical scenarios for the 2D 
THAFM.^'^-^ It is also possible that the paramagnetic- 



state ESR line is statically broadened, since neutron scat- 
tering in NiGa2S4^^'^^ and muon relaxation in both com- 
pounds suggest inhomogeneous spin structure above T* 
that could be reflected in the susceptibility. Measure- 
ments of the ESR linewidth field dependence in the para- 
magnetic state might clarify the situation. There is no 
evidence for a crossover above T* in the muon relaxation 
data for NiGa2S4, which agree quantitatively with Eq. (4) 
up to ~7 K above the ESR anomaly temperature. 



B. T < T*: quasistatic relaxation 

Quasistatic muon spin relaxation (the onset of nonzero 

(Bloc)) is observed in /iSR experiments below T* in both 
NiGa2S4 and FeGa2S4. An important issue is whether 
or not the quasistatic field itself is truly static or slowly 
fiuctuating. We discuss this in terms of two limiting sce- 
narios: (1) a static (Bioc) together with a fluctuating 
component (5Bioc(t), and (2) slow fluctuations of (Bjoc) 
as a whole. We designate these scenarios "conventional" 
and "unconventional" , respectively, for reasons that will 
become apparent, and consider the correlation time Tc as- 
sociated with (Bloc) • If quasistatic relaxation is observed 
Tc must at least be longer than l/\d (and possibly infl- 
nite), since otherwise motional narrowing would result in 
a single muon relaxation function rather than the two- 
term form of Eq. (l).^^ 

In the conventional scenario the dynamic muon relax- 
ation is due to transitions induced by (5Bioc(t) between 
muon spin Zeeman levels in a truly static (Bioc) (tc = 
oo). This is the normal situation in an ordered mag- 
netic phase, where S'Biodt) is due to thermally-excited 
magnons. These give rise to motionally-narrowed relax- 
ation, because the frequencies associated with spin-flip 
scattering of magnons are much higher than 'y^SBioc-'^^ 
In the unconventional scenario the miion relaxation is 
adiabatic,*5°'^'' and Tc « 1/^d ~ 10"^ and lO"'' s for 
NiGa2S4 and FeGa2S4, respectively, below T* [Fig. 3(a)]. 
This would be extremely slow on the time scale of the 3<i 
exchange interaction (~10~^^ s), and would correspond 
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to a highly unusual "spin-gel" phase as has been dis- 
cussed in the literature. ^"^'^'-^^ 

Unfortunately the /iSR data do not distinguish be- 
tween these scenarios, and the observed value of 1/Ad 
is therefore only a lower bound on Tc- Strong muon 
relaxation and NQR signal wipeout^'^ below T* could 
also occur in the conventional scenario/^ Nevertheless 
an extended critical regime or spin-gcl phase, where con- 
stituent spins arc strongly correlated but still fluctuat- 
ing on the microsecond time scale with zero long-time 
average, is compatible with the /itSR data for T < T*. 
The absence of a sharp anomaly in bulk properties at T* 
would be highly unusual in the conventional scenario but 
is perhaps not ruled out in the unconventional scenario; 
this feature may be evidence in favor of the latter. More 
work is needed to clarify the situation. 

As noted in Sec. I, the observation of the transition 
at T* in ^^Fe ME experiments is evidence against the 
conjecture^'^ that the muon electric charge perturbs its 
environment and somehow induces the transition ob- 
served in /iSR. In NiGa2S4 substitutional Fe^"*" and host 
Ni^+ ions are isovalent, and in Nio.99Feo.oiGa2S4 the ME 
data yield a value of T* close to that from /iSR in the un- 
doped end compound. Fe^"*" is of course itself the host 
in FeGa2S4. 



C. Comparison with other models 

We compare the /xSR results in FeGa2S4 and NiGa2S4 
with candidate models for 2D THAFM behavior other 
than CSS. 

Z2 vortex binding-unbinding transition. Kawamura 
and co-workers^^° and others'"^ have carried out exten- 
sive studies of Z2 vortices as topological defects in a 2D 
THAFM with nearest-neighbor interactions. Z2 vortices 
been studied in detail only in the classical limit, but pre- 
sumably they also exist for quantum spins. It was sug- 
gested that binding of thermally-excited Z2 vortices with 
decreasing temperature could lead to an unusual thermo- 
dynamic phase at a transition temperature T„, with long 
but finite spin correlation lengths and times below T^,.^'^" 

It has been conjectured*"^"'^® that the onset of qua- 
sistatic muon spin relaxation at T* might arise from this 
transition, so that T* = T^. The existence of a crit- 
ical regime between Tf and T* has been proposed for 
NiGa2S4 with Tf/T* - 0.4 (Refs. 18 and 29) and for 
NaCr02 with Tj/T* ~ 0.3.^^ From the present results in 
FeGa2S4 Tf/T* « 0.5, although we note again that the 
signature of Tf is very different in the susceptibilities of 
FeGa2S4 and NiGa2S4. The observed rapid muon spin 
relaxation below T* and the lack of a divergence of Ad at 
T* seem consistent with an important qualitative feature 
of the Z2 vortex picture, viz., finite spin correlation times 
at and below Ty. 

A number of questions remain concerning the /iSR re- 
sults and predictions of the Z2 vortex binding model: 



• In the model T-^ is slightly lower than the temper- 
ature of the lower specific heat peak, as is the case 
for T* in NiGa2S4.^'^ As previously noted,*^ how- 
ever, T* in FeGa2S4 lies between the two specific 
heat peaks, leading to doubt as to the applicability 
of the model to this compound. 

• The theory predicts only weak essential singulari- 
ties at Ty, where the fluctuation spectrum is domi- 
nated by the spin- wave contribution^' so that the 
CSS calculation should apply. This varies expo- 
nentially (i.e., rapidly) with temperature;, but a 
further abrupt increase of Tc on cooling through 
T* (Ref. 18) does not seem to be found in the Z2 
vortex picture; Tc merely increases (exponentially) 
with cooling through T^.^'^° In this case the maxi- 
mum rate is obtained for w^Tc « 1 (i.e., not neces- 
sarily at Ty),^^ and is of the order of > 200 fis~^, 
whereas the observed maximum rates in FeGa2S4 
and NiGa2S4 [Fig. 3(a)] are one to two orders of 
magnitude slower than this. Furthermore, such a 
slowing-down would not account for the mean-field- 
like temperature dependence oiijJij,{T) below T* ob- 
served in NiGa2S4,^^'^^'^^ since there would be no 
temperature dependence once cj^Tc ^ 1. Thus the 
data seem to rule out a smooth increase of Tc on 
cooling with no other eff'ects. 

• In the Z2 vortex model, spin dynamics above 

Ty involve both unbound vortices and conven- 
tional spin-wave excitations. A crossover from 
spin- wave-dominated to vortex-dominated spin dy- 
namics with increasing temperature is predicted 
at a temperature slightly above T^,.^'^° No such 
crossover or other sign of spin dynamics beyond 
the CSS prediction is observed in the /iSR data 
for NiGa2S4. As noted above a candidate for a 
crossover has been observed in ESR experiments 
on this compound, ^^'^^ but the applied field used 
in these experiments (20 T) was much greater than 
the value (^0.1 T) sufficient to change the symme- 
try of the Hamiltonian and suppress Z2 vortices.^ 

• An exponentially growing Z2 vortex density at 
high temperatures leads to a decrease in Tc.^'^° 
The resultant motional narrowing could explain 
the observed high-temperature dropoff of Xd/T^ 
in FeGa2S4 (Fig. 5) below the CSS prediction 
[Eq. (4)]. The curve in Fig. 5 is a fit to the FeGa2S4 
data of an expression of the form 

X, = A^^ (7) 

suggested by the Z2 vortex model, where Tsw and Ty 
are the spin correlation times associated with spin- 
wave and .Z2-vortex fiuctuations, respectively.^ 
Each correlation time incorporates its predicted 
temperature dependence, i.e., Tgw oc T^exp(To/T) 
and Ty (X exp{[Tz,/(T-T„)i"}, a « 0.5.^'^° The 
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fit is good (curve in Fig. 5), but the parame- 
ters are not well determined from the fit: = 
(1200 ± 2600) K and a = 0.5 ± 1.8. As noted in 
Sec. Ill A, however, the dropoff in FeGa2S4 may 
be due to other mechanisms, and cannot be un- 
ambiguously attributed to muon relaxation by Z2 
vortices. 

Thus there are points of agreement and points of dis- 
agreement between the Z2-vortex model and the experi- 
mental /zSR results. Additional work is necessary to re- 
solve these issues. In particular, a quantitative estimate 
of the expected muon relaxation rate due to unbound 
vortices above is necessary to determine whether the 
negative result in NiGa2S4 would be expected. 

Spin-nematic ground state, impurity spin dynamics. 
Treatments of the unconventional properties of NiGa2S4 
via "spin nematic" models with nonmagnetic quadrupo- 
lar ground states^^'^^'^'''"'^''''' are in agreement with the 
field-independent specific heat, but appear to be called 
into question by the observation of strong spin dynamics 
at low temperatures.^^ These theories have specifically 
considered only S = 1 systems,''^ and therefore do not 
explicitly address the similar behavior of the two com- 
pounds. However, the observation that the field indepen- 
dence and low-temperature behavior of the specific 
heat is preserved for integer-spin impurities in NiGa2S4 
but not for half-integer-spin impurities^^'^^ indicates the 
importance of integer spin for these properties. 

A recent theory by Takano and Tsunetsugu^^ con- 
cludes that bond disorder generates magnetic impurities 
in an antifcrroquadrupolar"'^"'^ spin-nematic state of a S=l 
2D THAFM with biquadratic interactions. Coupling to 
nonmagnetic excitations gives rise to an indirect long- 
range interaction between these impurities. "Freezing" of 
impurity spins occurs to a state with slowly-fluctuating 
spin moments with vanishing spin expectation value. The 
time scale of these fluctuations has not been reported. It 
is argued that slow fluctuations are also induced in the 
lattice spins, leading to a fluctuating field throughoiit the 
system for impurity concentrations of the order of 1%. 
Vortex unbinding strongly suppresses the impurity-spin 
interactions by preventing the definition of a consistent 
spin quadrupole director over a path around a vortex. 
This strong interplay between impurity spins and vortices 
leads to identification of T„ as the transition temperature 
for impurity-spin-moment freezing as well as vortex bind- 
ing. 

Like the .^2-vortex scenario, the spin-nematic impu- 
rity model is qualitatively consistent with a number of 
features of the /xSR data. The abruptness of the transi- 
tion, the possibility that the fluctuations are slow, and 
their suppression by relatively low magnetic fields coex- 
isting with a field- independent specific heat are all 
consistent with the /xSR results. Impurities give rise 
to magnetic inhomogeneity in the surrounding lattice, 
which might account for the observed sub-exponential 
muon relaxation. 



The nature of the impurity-spin fluctuations in this 
model is an important question. They appear to be quan- 
tum fluctuations, reflecting the vanishing of the expecta- 
tion value of the impurity spin in the ground state. One 
must be cautious in assuming that such ground-state fluc- 
tuations will result in probe-spin relaxation at a neigh- 
boring site, however, since the ground state is by defini- 
tion an eigenstate that does not fluctuate; thermal ex- 
citations are generally required for spin-probe relaxation 
in magnetic systems. 

As a simple example, consider two antiferromagnet- 
ically-coupled spins A and i?, where a probe spin is 
hyperfine-coupled only to spin A. One might expect the 
quantum fluctuations (oscillations) of A to induce tran- 
sitions between probe-spin states. It is straightforward 
to show, however, that the matrix elements of the spin 
operator Sa in the hyperfine coupling vanish in the sin- 
glet ground state, so that the probe spin is a constant of 
the motion even though the spin of A is not. It would 
be useful to clarify whether impurity-spin fluctuations in 
the spin-nematic impurity model are similarly ineffective 
at probe-spin relaxation or, alternatively, involve thermal 
excitation of degenerate or nearly-degenerate spin states. 
In the latter case persistent probe-spin relaxation at low 
temperatures would be expected. 

There are other questions concerning comparison of 
this theory with the experimental situation: 

• Consideration of a number of properties of NiGa2S4 
within the theory has not yet been reported. These 
include a quantitative estimate of the impurity- 
spin fluctuation rate, its temperature dependence, 
the origin and effect of spin freezing at Tf < T* , 
and the absence of persistent low-temperature re- 
laxation in NQR experiments, which is in marked 
contrast with its presence in the /iSR results. 

• Defect concentrations in NiGa2S4 have been found 
to be much lower than 1%,^^'^^ and /^SR experi- 
ments have been carried out on a number of sam- 
ples of NiGa2S4 from several laboratories with sim- 
ilar results. ^^'^^"^^ This suggests that impurities 
might not play a major role unless the concentra- 
tion dependence is expected to be weak. On the 
other hand, the observation of two resonances in 
gallium NQR experiments on NiGa2S4^" suggests 
structural inhomogeneity at some level; stacking 
faults have been suggested. 

• The identification of impurity-spin freezing with 
vortex binding below T„ is based on the pres- 
ence of unbound quadrupolar vortices above this 
temperature.^^ As discussed above, there is no evi- 
dence for unbound vortices above T* from /iSR 
in NiGa2S4, and it has been argued^ that Z2 vor- 
tices based on noncollincar AF order are the most 
likely species in this compound. Muon spin re- 
laxation by nonmagnetic quadrupolar vortex exci- 
tations would, however, be weak and likely to be 
masked by relaxation due to spin fluctuations. 
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We conclude that the spin-nematic model with impurity- 
spin freezing is also a promising candidate mechanism, 
but that more work is required to determine the extent to 
which it is supported by fiSR and other experiments. One 
important aspect is the need for evidence for or against 
the substantial biquadratic spin interaction that is nec- 
essary for a nematic ground state. 

Broken C3 symmetry. Models with nearest- and third- 
nearest-neighbor exchange couplings, with^'"^ or without 
biquadratic coupling, (exhibit ground states with broken 
C3 rotational symmetry. Growth of quadrupolar short- 
range order at high temperatures, proposed^^ to account 
for the high-temperature specific heat peak, is not easily 
tested by the /zSR data, since as noted above quadrupole 
excitations would be hard to detect with a magnetic 
probe such as the muon spin. The magnetic correla- 
tions that generate T* are associated with the lower spe- 
cific heat peak,^^ however, whereas in FeGa2S4 T* is well 
above this peak.^^ 



IV. CONCLUSIONS 

/iSR studies of FeGa2S4 indicate drastic slowing and 
possible freezing of magnetic fluctuations below an un- 
conventional transition at T* « 31 K, which is twice the 
spin glass-like freezing temperature T/ w 16 K from mag- 
netic susceptibility measurements. Muon spin relaxation 
rates above T* in both FeGa2S4 and NiGa2S4 are in very 
good quantitative agreement with the CSS result^^ for 
spin-lattice relaxation in a 2D quantum antiferromagnet. 
An extended regime of strong spin fluctuations is ob- 
served for T < T*, and the spin dynamics scale between 
the compounds from ^Tj to above T*. 

Although there are differences in the magnetic prop- 
erties of FeGa2S4 and NiGa2S4 (most notably the very 
difii'erent anomalies at T/), similarities in the /uSR data 
include the 2D critical spin dynamics above T*, the ab- 
sence of a divergence of at T* , and scaling between the 
two compounds in the neighborhood of T*. These fea- 
tures are qualitatively consistent with theories of phase 
transitions driven by two distinct types of defects; (1) Z2- 
vortex topological excitations, and (2) magnetic im- 
purities in a nonmagnetic spin-nematic ground state. 
There are questions concerning reconciliation of either 
model with other aspects of the data, however. The few 
differences in the ^SR behavior of the two compounds 
include a possible anomaly in Ad(T) at Tf in NiGa2S4 
[Fig. 3(a), but see the associated discussion] and the 
dropoff of Xd{T)/T^ at high temperatures in FeGa2S4 
(Fig. 5). 

It should be noted that in general an exponential 

temperature dependence is expected from any large- 
scale excitation with constant energy. As discussed in 
Sec. Ill A, the strongest evidence that in FeGa2S4 and 
NiGa2S4 muon relaxation is dominated by the specific 
spin-wave excitations treated by CSS is the remarkable 



agreement of the /iSR values of the exchange constants 
with those from other experiments (Table I). 

The /iSR experiments do not discriminate between 
slow fluctuations and truly static freezing below T*. Nev- 
ertheless, a sharp anomaly is virtually always seen in 
some bulk property at a magnetic transition. In FeGa2S4 
the absence of an anomaly at T* in the susceptibility and 
at Tf in the muon spin relaxation appears to rule out a 
transition to a fully spin-frozen state at cither of these 
temperatures. This and the coexistence of strong and 
strongly field-dependent magnetic spin dynamics with a 
field-independent specific heat at low temperatures are 
probably the two most remarkable results of this study. 
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Appendix: Spin-lattice relaxation in the 
renormalized-classical and quantum-critical regimes 
of a 2D frustrated antiferromagnet 

We consider the CSS calculation of the spin-lattice re- 
laxation rate l/7i (Ref. 49) in more detail, and com- 
pare their results in the renormalized-classical (RC) and 
quantum-critical (QC) regimes of a 2D frustrated antifer- 
romagnet that orders at T = 0. We follow the equation 
numbering and notation of CSS. 

In the RC region 



1 



(RC) 



Ao 

h 



) , keT «; 27rp« 
Ps J 



(4.29) 



where ^0 is the hyperfine coupling constant, A^o is the 
condensate magnitude (renormalized by quantum fluctu- 
ations but presumably of order unity), ^ is the correlation 
length, c is the spin-wave velocity, and Ps is the spin- 
wave stiffness constant. The prefactor required to make 
Eq. (4.29) an equality is of order unity.^® The correlation 
length in the RC region is given by 



^^1^ he ({N-l)kBT 



2^ knT 



4nps 



l/2(iV-l) 



exp 



where A'' = 2 is the dimensionality and 



{N-l)kBT\ ' 
(4.7) 



l/2(Ar-l) 



X r(l -M/2(iV-l)) (4.E 



11 



xr(3/2)=(ie)'/'^ = 0.5166. 

(A.1) 



The CSS result for 1/T^^''~^'' is proportional to the ra- 
tio ^/c. In Sec. Ill A we note that the discrepancy be- 
tween the fit value of NqAq/Ti in NiGa2S4 and the ex- 
pected value ^a;p(0) might be related to uncertainty 
in this quantity. Neutron scattering experiments in 
NiGa2S43i obtained ^ = 25 A at 10 K and spm-wave 
velocity he = 29 meV A at 1.5 K, yielding ^/c = 
5.7 X 10"-^^ s. Using the value of ps for NiGa2S4 from 
Table I, the calculated CSS value of ^/c from Eq. (4.7) at 
T* = 9 K is 1.04 X 10~^ s, more than five orders of mag- 
nitude larger than the neutron scattering value. Assum- 
ing NoAo/h Ri a;p(0), the relaxation data and Eq. (4.29) 
yield ^/c « 2.5 x 10~^° s, more than two orders of mag- 
nitude smaller than that from Eq. (4.7) but still much 
larger than the neutron scattering value. Clearly this 
parameter is not well understood at present. 
In the QC region the CSS result for 1/Ti is 



1 



rpiQC) \ h J Ps y 4:Trps 

(5.15) 

where Z = {^5 - 1)/8N = 0.07725 for N = 2, and 
1] = 1 + 32/37r^A'' is a scaling exponent = 1.5404 for 



n 



NksT 



N = 2. 

To compare the results in the RC and QC regions we 
define the dimensionless quantities 



1 Ps I f h 



1 



, t = kBT/2np, , (A.2) 



so that 



= 18.074 exp(2/t) , t <^ 1 , (A.3) 



and 



,(QC) 



= Zt'> 



= 0.07725 i>l. 



(A.4) 



In Fig. 6 l/tit^ is plotted versus 1/t for each region all 
the way to 1/t = 1 [i.e., beyond the validity of Eqs. (A.3) 
and (A.4)]. It can be seen that at 1/i = 1 the extrap- 
olated value of 1/tit^ is much larger for the RC regime 
than for the QC regime. Thus in an exact solution 1/ti 
would be expected to drop below the asymptotic RC re- 
sult [Eq. (A.3)] as the QC region is approached, although 
the caveat concerning the prefactor in Eq. (4.29) should 
be noted. 
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